Quantum teleportation of optical images with frequency conversion 
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We describe a new version of continuous variables quantum holographic teleportation of optical 
images. Unlike the previously proposed P, Q scheme, it is based on the continuous variables 
quantum entanglement between the light fields of different frequencies and allows for the wavelength 
conversion between the original and the teleported images. The frequency tunable holographic 
teleportation protocol can be used as a part of light-matter interface in quantum parallel information 
processing and parallel quantum memory. 
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Continuous variables quantum teleportation allows to 
transfer an arbitrary quantum state of the electromag- 
netic field between two spatially separated systems via 
an exchange of classical information in combination with 
quantum entanglement shared by these systems. Initially 
the protocol of continuous variable teleportation was pro- 
posed theoretically 0,0] and realized experimentally[a,0| 
for a single spatial mode of the electromagnetic field. 
Recently this protocol was generalized theoretically to 
spatially multimode electromagnetic fields 0, Q. The 
important feature of spatially multimode teleportation 
scheme is that it allows for simultaneous parallel telepor- 
tation of optical images (still or time- varying) containing 
large number of elements or pixels. The proposed scheme 
was called holographic teleportation because it resembles 
the conventional holography with an important differ- 
ence that the reconstructed image is a quantum copy of 
the original one with fidelity that can be made close to 
unity. It is impossible to achieve high fidelity in holo- 
graphic teleportation without sharing multimode quan- 
tum entanglement by the original and the target systems. 

One of possible applications of quantum teleportation 
is the light-matter interface for quantum memory^] for 
light which allows to record and story a quantum state 
of light on that of an atomic ensemble. The generalized 
spatially multimode teleportation scheme opens new per- 
spectives for creation of a parallel quantum memory for 
parallel processing of quantum information. For an effi- 
cient interaction between the light wave and the atomic 
medium it is desirable to have a possibility for tuning 
the optical frequency of light to that of the atomic tran- 
sition without changing the quantum state of the elec- 
tromagnetic field. Let us mention a proposal of quantum 
frequency conversion!^ that was the first demonstration 
of quantum state transfer into different frequency. In 
recent experiment quantum transfer of qubits was suc- 
cessfully demonstrated between photons of wavelength 
1310 nm and 710 nm. 

In this letter we describe a new version of holographic 
teleportation of optical images with frequency conversion. 



This means that an input image, carried by a light wave 
with optical frequency toi, is teleported to the output 
image at different optical frequency W2 with preservation 
of its original quantum state. Similar to the frequency- 
preserving quantum holographic teleportation, discussed 
in Refs. 0j 0i to obtain high fidelity in the teleported 
image one needs to create spatially multimode entangle- 
ment between the original and the target systems. In 
our case the source of spatially multimode entanglement 
between the input image at frequency uji and the tele- 
ported image at frequency uj2 is a type-I traveling-wave 
non-degenerate optical parametric amplifier (OPA). 

The optical scheme of holographic teleportation with 
frequency conversion if shown in Fig. 1. An input im- 
age which is to be teleported from Alice to Bob, is de- 
scribed by a slowly- varying field operator A™ (p, t) , where 
p— {x,y) is a two-dimensional transverse coordinate. An 
input electromagnetic wave carrying this image is splitted 
into two secondary waves by a 50/50 beam splitter BSi, 
and two quadrature components of these secondary waves 
are homodyne detected by two difference detectors using 
two local oscillators LOa; and LOy, two 50/50 beam split- 
ters BS2 and BS3, and four efficient CCD cameras with 
appropriate spatial resolution. The local oscillator waves 
have the same frequency ui as the carrier frequency of 
the input image. The slow-varying difference photocur- 
rent densities from these CCD cameras, containing infor- 
mation about spatio-temporal quantum fluctuations of 
the input image, are transmitted from Alice to Bob and 
are used for preparation of the output field A2{p, t). The 
multi-channel modulators M^; and My perform spatio- 
temporal modulation of a coherent wave with different 
carrier frequency lu2. This part of the scheme can be 
viewed as a classic non-stationary holography, where the 
reconstructing wave has different wavelength. 

A key ingredient of our teleportation scheme is a pair of 
spatially-multimode EPR fields En{p,t), n — 1,2. Since 
an entanglement between different carrier frequencies is 
needed, the fields are created by a type-I traveling-wave 
non-degenerate OPA. The dichroic mirror Mi reflects the 
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FIG. 2: Uncertainty region for the slow field amplitude from 
Eq. 4 (a), the frequency dispersion of squeezing ellipses for 
non-degenerate traveling- wave OPA, where r2(0, 0) = 3, - 
arb. units (b). Here X, Y - field quadrature components. 



FIG. 1: Optical scheme of holographic teleportation with fre- 
quency conversion. 

wave at frequency uji and transmits that at frequency 
102- In the case of perfect reflectivity there are no fre- 
quency matched vacuum fluctuations entering from the 
open ports of Mi into the corresponding fields. 

Let us introduce slowly-varying spatio-temporal an- 
nihilation and creation operators En{p,t) and E^^{p,t), 
n = 1,2 of the electromagnetic waves with central fre- 
quencies tui and UJ2 at the output of the non-degenerate 
OPA. The frequencies cui and uj2 obey the condition of 
energy conservation, uji + uj2 = ojp, where tOp is the fre- 
quency of the pump wave. The field operators are nor- 
malized so that {E}^{p, t)En{p, t)) gives the mean value of 
the irradiance, expressed in photons per cm^ per second. 

The transformation of the input fields An (p, t) of the 
non-degenerate OPA in the vacuum state into the out- 
put fields En (/O, t) in the broadband multimode squeezed 
state is described in terms of the Fourier components of 
these operators in frequency and spatial-frequency do- 
main, E{p,t) e(g, fi). The squeezing transformation 
performed by a non-degenerate OPA, can be written as 
follows: 

eniq,n) - Un{q,n)an{q,n) + Vn{q,n)ali-q,-n), (1) 

where n = 1,2, n 7^ n' , the coefficients Un{q,^) and 
Vn{q, i^) depend on the amplitude of the pump field, non- 
linear susceptibility and the phase-matching condition. 
The explicit form of UjJq, n) and Vn{q, ^) can be found, 
for example, in Ref. [llj . The EPR correlations between 
the fields Ei (p, t) and E2 {p, t) are determined by two pa- 
rameters, namely, the orientation angle ipniq-,^) of the 
major axis of the squeezing ellipse |lO| and the degree of 
squeezing r„((f, 17), 

Mq,^) ^lsiTg{Uniq,n)Vn'i-q,-n)}, (2) 



exp [±r„(<r, n)] = \Un{q, n)\ ± \Vn'{q, fl)]. (3) 

It should be noted that in Fig. 1 we have chosen the 
input field at the frequency uji and the output field at 
the frequency 0^2, that corresponds to n = 1 and n' — 2 
in Eqs. (O, ©, but with the same source of entangled 
beams the teleportation with frequency conversion can 
be performed from 102 to cui as well. 

It is instructive for better understanding of the scheme 
to give a simple physical explanation for the continu- 
ous variables quantum entanglement arising from squeez- 
ing in non-degenerate parametric down-conversion. For 
simplicity we shall omit the spatial dependence of the 
fields. Let us introduce the sum of two slow amplitudes, 
£{t) ^ Ei{t) + E2(t) and consider a contribution to £{t) 
from a pair of Fourier amplitudes, 

£{t) ~ ei(-r!)e'^^* H- e2[n)e~'^\ (4) 

One can show [fj], that the discussed above squeezing 
ellipse for given f2 represents exactly the quantum uncer- 
tainty region for this contribution. For perfect squeezing 
this automatically implies a perfect phase matching be- 
yond the classical limit between the slow amplitudes in 
the right side of as illustrated in Fig. 2a. Let us 
remind that these amplitudes were defined with respect 
to the different carrier frequencies loi and UJ2- Hence 
the corresponding beams can be split by a dichroic mir- 
ror, preserving their quantum entanglement. The first 
slow amplitude is imprinted into the photocurrents by 
the homodyne detection and into the reconstructing field 
at different frequency UJ2 by the modulation. The sec- 
ond one ab initio is related to this frequency. A proper 
matching of the relative amplitudes and phases in quan- 
tum and classical channels leads to the quantum noise 
cancellation and effective teleportation. 

Next important point is that the quantum modulation 
of the output field at a given frequency is now due to two 
independent pairs of entangled amplitudes, the first pair 
is shown in (0}, another one corresponds to — s- — fi. 
Therefore, in order to achieve optimal results one has to 
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take care of a proper phase matching of two indepen- 
dent squeezing elhpses (at V, and —ft). The frequency 
dependence of squeezing within the frequency range of 
effective non-degenerate parametric down-conversion is 
graphicaUy shown in Fig. 2b. 

Apart from these new features, introduced by the fre- 
quency non-degenerate nature of the process, the fields 
evolution in the scheme is very similar to that described 
in Refs. The input field 4'"(p, t) is mixed with one 

EPR beam Ei{p,t) at the 50/50 beam splitter BSi. The 
secondary waves and By, 

Bx,y {p, = ^ ( ± t) + El {p, t)) , (5) 

with the +(— ) sign corresponding to x{y) channel, are 
photodetected by means of two balanced homodyne de- 
tectors. The difference photocurrents Ix{p, t) and Iy{p, t) 
collected from individual pixels of these matrices, carry 
the information about the spatio-temporal quantum fluc- 
tuations of the quadrature components of Bx{p,t) and 
By{p, t), phase matched with LO^^ and LOy. These pho- 
tocurrents are sent from Alice to Bob via two multichan- 
nel classical communication lines and are used by Bob for 
spatio-temporal modulation of an external coherent wave 
with the frequency ^2, by means of two multichannel 
modulators and My. The teleported field A'^''\p, t) is 
created by mixing of this modulated wave with a second 
EPR wave E2 {p, t) at the mirror M with high reflectiv- 
ity. By choosing appropriately the mirror transmission 
and the modulation depth, the teleported field (p, t) 
can be obtained as 

AT'{p,t) = A^,^ip,t)+F{p,t), (6) 

where F{p,t) is an operator describing the noise added 
to the teleported field, 

F{p,t)^E2{p,t) + El{p,t). (7) 

It can be demonstrated that the added noise can be con- 
sidered as classical, thus preserving the commutation re- 
lations for the output field. 

An ideal teleportation from oji to UJ2 would correspond 
to the situation when the output field operator at dif- 
ferent spatial points p and at different time moments t 
is an exact copy of the input field operator, A™* (/O, t) = 
A™{p, t). However, as explained in Refs. Qilll, this would 
require an infinite energy of the EPR beams and therefore 
could never be realized in practice. In real experimental 
situation the teleportation process will never take place 
"point-to-point" , but on average within some finite spa- 
tial area and within some finite time interval. In order to 
quantitatively describe the performance of our telepor- 
tation scheme we introduce a coarse-grained description 
of the input and output fields. Precisely, we consider 
the averaged field operator over a square pixel Sj of area 



S' = and over a time interval Ti of duration T: 

ir(j,*) = i/ dp I dtAr'{p,th (8) 

Vol Js- JTi 

with analogous definitions for the input field. The av- 
eraged field operators obey standard commutation rela- 
tions of discrete field oscillators. The quadrature compo- 
nents of the averaged field operators at the output are 

x;^\j,i) = ir(j,*)e-'^ + /i.c., 

= -^Ar\J,^)e''^ + h.c., (9) 

and similarly for the input field operators. 

To characterize the noise added in the teleportation 
process we compare the correlation functions of the input 
and output quadrature components defined in Eq. 
As follows from Eq. ©, the relation between these cor- 
relation functions is 

(5X°"'(j,*)JX™*(j',i'))= (10) 

where C{j, is the covariance matrix of the added 

noise. The explicit expression of this covariance matrix 
is 

CU,f;i,i') = 2 J dqj dnBA{q)BT{n) x (11) 

cos [q - {pj - pj,) - n{u - ti,)]G{q,n), 

where pj = {xj, yj} is the center of the j-th pixel, and ti is 
the center of the i-th time interval. The functions B/\{q) 
and Bt{^) arise from the coarse-graining operation, and 
for a square pixel of size A they read, 

BT(fi) = §sinc^ (^) . (12) 

The Green function G{q, fl) expressed in terms of the 
orientation angle ■02(9, ^^) and the degree of squeezing 
r2 ((f, il) from Eq. ^ looks like 

Giq, n) = e2''^('"'") cos^ ij2iq, r!) + e-2'-2(9-,n) ^^^^2 ^^^^^ 

(13) 

As follows from Eqs. (|11|I - H13I) . the covariance matrix 
C(j, is independent of the phase (p in Eq. 0. 

Thus, the added noise is the same for any quadrature. 

In the absence of the EPR correlations, i. e. when 
r2((f, n) = 0, we obtain the classical limit of teleporta- 
tion with the covariance matrix Cci{j,j';i,i') = 25jj>diii. 
In this limit two units of vacuum noise are added at each 
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FIG. 3: Diagonal elements C{j,j;i,i) of the noise covariance 
matrix as a function of the relative pixel size D = A/lc for 
three different observation times T equal to lOTc, Tc, and 
O.lTc, here r2(0,0) = 3. 



pixel exactly as for a single-mode teleportation and for 
the teleportation of images without frequency conversion. 

Choosing the phase matching condition in the OPA 
such that ^/'2(0, 0) = n/2, we obtain reduction of the 
quantum noise below the standard quantum level within 
some bandwidth of spatial frequencies q and temporal 
frequencies fl determined by the phase-matching condi- 
tions in the crystal. When the pixel size A and the time 
window T are much larger than the characteristic coher- 
ence length Ic and the coherence time Tc of the OPA, we 
obtain 



lim 

+ 00 .T- 



Cc\{j,j';i,i') 



exp[-2r(0,0)]. (14) 



In Fig. 3 we illustrate the role of the pixel size A and of 
the integration time T in the teleportation process with 
frequency conversion. This figure shows the diagonal el- 
ements C{j,j]i,i) of the covariance matrix as a function 
of the relative pixel size D — A/lc for three different 
observation times T equal to lOTc, Tc, and Q.lTc- The 
coherence time Tc for a frequency non-degenerate OPA 
is typically estimated as the time delay at the crystal 
length / between two wave packets centered at the fre- 
quencies uji and 102, Tc l\i/vi — 1/^2 1 : arising due to 
the difference of the group velocities Vn- Since a typ- 
ical frequency spectrum of parametric down-conversion 
is fairly broad, a wave-packet spread due to the group 
velocity dispersion can also have an effect on squeezing 
and entanglement in our teleportation scheme, as we shall 
discuss elsewhere [l^. 

The frequency dependence of squeezing orientation an- 
gle, arising in a parametric crystal and illustrated in 
Fig. 2b, is a pure phase effect and can be compensated by 
propagation in a linear medium with a properly chosen 
frequency dependence of the refraction index. It is worth 
noting that in many interference experiments with twin 
photons, a need for similar compensation has been real- 
ized some time ago. Such a compensation can be applied 
to the plots shown in Fig. 3 and will be illustrated in the 
forthcoming publication [l^ . 

The coherence time for a non-degenerate OPA is usu- 
ally much larger than for a degenerate one. Therefore, 



the assumption T ^ Tc used in Refs. 0, is not neces- 
sarily true anymore. The role of the observation time T 
is as crucial in the teleportation process as the role of the 
pixel size A. Indeed, as follows from Fig. 3, the diagonal 
elements of the covariance matrix which characterize the 
added noise power, decrease with growing pixel size A 
from the classical limit C{j,j;i,i) = 2 to the EPR limit 
= 2exp[— 2r(0,0)] when the observation time 
is large compared with the coherence time Tc- However, 
when the observation time T becomes comparable with 
the coherence time Tc or smaller, the diagonal elements 
j; i, i) never reach the EPR limit but remain close to 
the classical limit even for large pixel size. 

This results allow to estimate an effective number of 
spatio-temporal degrees of freedom for an input non- 
stationary image that can be teleported in our scheme. 

In conclusion we have proposed a new scheme of quan- 
tum holographic teleportation of optical images from one 
optical frequency to another. We have presented pre- 
liminary analysis of the performance of our scheme and 
have revealed the possibilities to achieve high quality tele- 
portation. More detailed calculations including evalua- 
tion of multi-pixel fidelity for this scheme are in progress 
and will be published elsewhere. This work was sup- 
ported by the Network QUANTIM (IST-2000-26019) of 
the European Union and by the INTAS under Project 
No. 2001-2097. The research was performed within the 
framework of GDRE "Lasers et techniques optiques de 
I'information" . 



[1] 

[2] 

[3] 
[4] 

[5] 
[6] 

[7] 

[8] 

[9] 

[10] 
[11] 

[12] 



I. V. Sokolov, M. I. Kolobov, A. Gatti, and L. A. Lugiato, 

Optics Communications, 193, 175 (2001). 

A. Gatti, I. V. Sokolov, M. I. Kolobov, and L. A. Lugiato, 

Eur. Phys. J. D 30, 123 (2004). 

L. Vaidman, Phys. Rev. A 49, 1473 (1994). 

S. L. Braunstein and H. J. Kimble, Phys. Rev. Lett. 80, 

869 (1998). 

A. Furusawa, J. L. Sorensen, S. L. Braunstein, 
C. A. Fuchs, H. J. Kimble, and E. S. Polzik, Science, 
282, 706 (1998). 

W. P. Bowen, N. Treps, B. C. Buchler, R. Schnabel, 
T. C. Ralph, H.-A. Bachor, T. Symul, and P. K. Lam, 
Phys. Rev. A 67, 032302 (2003). 

B. Julsgaard, J. Sherson, J. Fiurasek, J. I. Cirac, and 
E. S. Polzik, Nature 432, 482 (2004). 

J. Huang and P. Kumar, Phys. Rev. Lett. 68, 2153 
(1992). 

S. Tanzilh, W. TiUel, M. Haider, O. Alibart, P. Baldi, 
N. Gisin, and H. Zbinden, Nature 437, 116 (2005). 
M. I. Kolobov, Rev. Mod. Phys. 71, 1539 (1999). 
E. Brambilla, A. Gatti, M. Bache, and L. A. Lugiato, 
Phys. Rev. A 69, 023802 (2004). 

L. V. Magdenko, M. I. Kolobov, and I. V. Sokolov (in 
preparation) . 



